Coupled magnetic and ferroelectric domains in multiferroic Ni3V208 
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Electric control of multiferroic domains is demonstrated through polarized magnetic neutron 
diffraction. Cooling to the cycloidal multiferroic phase of Ni3V208 in an electric field E causes the 
incommensurate Bragg reflections to become neutron spin polarizing, the sense of neutron polar- 
ization reversing with E. Quantitative analysis indicates the E-treated sample has handedness that 
can be reversed by E. We further show close association between cycloidal and ferroelectric do- 
mains through E-driven spin and electric polarization hysteresis. We suggest that definite cycloidal 
handedness is achieved through magneto-elastically induced Dzyaloshinskii-Moriya interactions. 

PACS numbers: 75.25.+Z, 75.60.-d, 75.80.+q, 77.80.-e 



Materials that are both ferroelectric and magnetic are 
classified as multiferroics. In some multiferroics where 
the ferroelectric and magnetic phases coexist, spin and 
charge are strongly coupled, leading to the possibility of 
controlling magnetic properties through an electric field 
E. Such a nonlinear magneto-electric response is of fun- 
damental interest and holds the potential for applications 
that include sensing, spintronics, and microwave commu- 
nication [1]. Recent studies have shown that an external 
E applied to multiferroics with non-collinear spin struc- 
tures, such as TbMnOs and LiCu202, favors a particular 
handedness of the magnetic order [ij, 0] • Other studies 
have shown the E control of domain population related 
to equivalent magnetic propagation vectors [4]. Here we 
examine the suppression and promotion of cycloidal mag- 
netic structures in Ni3V2 08 (NVO) by an applied E. Our 
quantitative analysis of the polarized magnetic diffrac- 
tion cross-section and hysteresis curve for this multifer- 
roic material indicates that a clockwise cycloidal single 
crystal can be generated and is stabilized by magneto- 
elastically induced Dzyaloshinskii-Moriya interactions. 

NVO is an insulating magnet with spin-1 Ni^+ ions ar- 
ranged in a buckled kagome-staircase geometry [5] . The 
spins occupy two distinct crystallographic sites denoted 
cross-tie and spine [See Fig. [T](a)]. Competing nearest 
and next-nearest neighbor interactions along the spines 
yield a complex magnetic phase diagram [6]. Magnetic 
inversion symmetry breaking was inferred in the so-called 
low-temperature incommensurate (LTI) phase, where un- 
polarized neutron diffraction data indicate a magnetic 
cycloidal structure with spins in the a-b plane and py- 
rocurrent measurements find concomitant electric polar- 
ization along the b axis. A Landau mean field theory 
was previously devised to account for this multiferroic 
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FIG. 1: (Color online) (a) NVO crystal sublattice showing 
Ni^^ spine (red) and cross-tie (blue) sites, (b) Counter- 
clockwise (top) and clockwise (bottom) spin cycloids prop- 
agating along the a axis. The (green) vertical arrow indicates 
the direction of P. 



behavior [sl, 0, Q • The free-energy expansion is 

F = a(T- ThWh + KT - Tl)cjI + 0{a^) 



V. 



(1) 



Here, a and h are constants, T is temperature, gh and 
cfl are the magnetic order parameters in the high-T in- 
commensurate and LTI phases, respectively, xe is the 
electric susceptibility, and P is the electric polarization. 
The last term is the lowest order (trilinear) symmetry- 
allowed multiferroic interaction, which in the LTI phase 
is given by Vlti = CL^o-HCfLP-f- Minimizing F with 
respect to P, one finds that P5 = b • P varies with T 
in proportion to the product of the two magnetic order 
parameters (P5 oc a^XEO'LO'H)^ as observed experimen- 
tally [9]. The theory also suggests that ferroelectric and 
magnetic domains are coupled in NVO. Here we examine 
this hypothesis by probing the magnetic and ferroelectric 
response to an applied E in the multiferroic phase. 

NVO crystals were grown from a BaO-V205 fiux 0. 
The buckled kagome layers span the a-c crystallographic 
plane and form the largest crystalline surfaces. A 0.58 
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FIG. 2: (Color online) Full polarized magnetic diffraction in 
an electric field under opposite spin-flip scattering conditions. 
Statistical uncertainties represent one standard deviation, (a) 
Cooling in a field E= +300 kV/m from 11 K to 5 K (E- 
cooling) favors cycloidal domains that predominantly diffract 
neutrons polarized antiparallel to Q. (b) Reversing E yields 
the opposite polarized intensity asymmetry. 



g, 120 mm^ crystal was selected for this experiment. A 
parallel-plate capacitor was formed by evaporating a 5 
nm Cr/40 nm Au layer on each large face of the crys- 
tal, hence normal to the ferroelectric axis. Au wires were 
attached to each side of the sample using silver epoxy 
paste. Polarized neutron diffraction measurements were 
carried out on BT-7 at NIST. A 14.7 meV neutron beam 
was polarized and analyzed by ^He neutron spin filters 
flQ ]. Helmholtz coils were used to generate a guide field 
at the sample position, thus defining the neutron spin 
quantization axis. In the vertical field (VF) configura- 
tion the field strength was 0.4 mT, oriented normal to 
the scattering plane to within 0.07 rad. In the horizon- 
tal field (HF) configuration, the field strength was 0.8 
mT, oriented within 0.1 rad of the horizontal plane and 
parallel to wave vector transfer Q = ki — kf to within 
0.2 rad. Mezei neutron spin- flippers were mounted be- 
fore and after the sample, providing access to a total of 
eight configurations for the incoming and outgoing neu- 
tron spin. The nomenclature used is as follows: (+) refers 
to a flipper off and (— ) refers to a flipper on. With both 
flippers off the neutron spin nominally points up for VF 
and parallel to Q for the HF configuration. 

HF spin-polarized diffraction data are shown in Fig. 
[21 There is a strong asymmetry in the intensity be- 
tween the (+/— ) and (— /+) spin- flipper configurations 
which reverses with E. This demonstrates in a quali- 
tative fashion the E-driven suppression and promotion 
of cycloidal magnetic domains [11]. Examining the data 
more carefully, we note that a finite peak remains under 
the (+/—)[(—/+)] spin-flip scattering condition, even af- 
ter cooling in a +300 [—300] kV/m field. Similar ef- 



TABLE I: Spin components on Ni^^ spine (s) and cross-tie 
(c) sites describing CW and CCW cycloids. The inversion 
operator T converts a CW cycloid into a CCW as follows: 
for a = a or c, X[mJ^] = T(^?i)* ^^le b-component, 

^\r^s%\ — i(^si)* (upper sign for i — V\ and lower sign for 
i — V^. For all cross-tie sites T\mc%\ — (vfici)* . There are 
six additional atoms in the conventional unit cell, obtained 
by translating di by (| , | , 0) . Asterisks denote complex con- 
jugation. 
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(i,-0.13, \) md, =(mJi,m^i,m^i) + (mj4,mS4,m^4) (nidg)* 



(i,0.13, I) 
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> =(Kl'^sl'^sl) + (^?4'^s4'^s4) (md4)* 

^ = (mJi,m^i,m^i) + (mj4, 771^4, m^4) (mdj* 

L = (^sl'^sl'^sl) + (^J4'^s4'K4) (md2)* 

md5 =(m«,,0,0) + (0,m^4,m-4) (nidj* 
mde = (m^i,0,0) + (0,m^4,m^4) (md^)* 



fects have been seen in [l, [sl, but have not been fully 
accounted for. Note that, because the strong diffraction 
cross-section is suppressed by both the incident and the 
final beam ^He neutron spin filters, the residual inten- 
sity cannot be accounted for by finite beam polarization. 
Quantitative analysis of the polarized diffraction cross- 
section is needed to account for this effect. 

Within a single domain, the magnetic structure of 
NVO in the LTI phase can be described as follows: 



5Rdi = nidie 



iqm-(R+di) 



iqm-(R+di) 



(2) 



Here, qm is the magnetic propagation vector, R is a vec- 
tor from the origin to the unit cell, are position vectors 
for Ni^+ ions within the conventional unit cell, and 
transform according to irreducible representations of the 
magnetic space group and specify the time-averaged mag- 
netic moments on Ni^+ sites. In the LTI phase, where 
electric polarization is present, the spin structure was 
previously described by the Fi and r4 irreducible repre- 
sentations with best fit basis vectors for spine and cross- 
tie sites m^i and iric^, i = 1,4 [7]. The resulting spin 
structure is a clockwise (CW) cycloid, progressing along 
a [see Fig. [T] (b)]. Spatial inversion is a symmetry op- 
eration of the paramagnetic phase that converts a CW 
cycloid into a counterclockwise (CCW) cycloid. The set 
of complex nidi and CCW cycloids are listed in 

Table [H We expect that domains in NVO are associated 
with these symmetry-related structures. 

To quantitatively characterize the cycloidal domains, 
we employed Blume's equations for elastic scattering of 
polarized neutrons [11]. The polarization-resolved differ- 
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FIG. 3: (Color online) Comparison between measured and 
calculated intensities for (1.27,1,1) and (1.27,3,1) peaks. 
Columns represent calculated intensities for single-domain 
(SD) and multi-domain (MD) cycloidal structures. (Red) cir- 
cles denote the average measured intensity and vertical bars 
show lowest and highest values measured for each particular 
condition. The categories include spin-flip (SF) and non-spin- 
flip (NSF) modes under HF and VF. + and — refer to the sign 
of Pi • E. Polarization efficiency corrections were applied in 
this comparison [T^]. 
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FIG. 4: (Color online) Evolution of magnetic and ferroelec- 
tric domains after E-cooling the sample. Ferroelectric data 
was obtained by integrating the pyroelectric current while E- 
heating to the paramagnetic phase at a rate of 1 K/min. (In- 
sert) Gold-plated NVO single crystal. 



ential diffraction cross- section is: 

+ n(T)[Pi-^±(T)]-Pi|^^(T)P 

Ti[^I(T) X^^(T)]}5(Q±q^-T). (3) 

Here, 7 =1.913, tq is the classical electron radius, N is 
the number of magnetic unit cells with volume vq^ i^Ni(T) 
is the form factor for Ni^+, r is the nuclear reciprocal 
lattice vector hosting the magnetic satellite, Pi and Pf 
are the initial and final neutron spin polarization, respec- 
tively, J-±{t) = Q X {J-{t) X Q), and the magnetic vector 
structure factor is given by 

Furthermore, the overall Pi-dependent diffraction cross- 
section is: 

all vq 

±zPi-[^l(T) x^^(T)]}(^(Q±q„,-r). (5) 

Dividing Eq. (|3]) by ([5]) provides the final polarization 
unit vector Pf. The calculated intensities for a single- 
domain (SD) cycloidal structure (CW or CCW) are read- 
ily obtained using Table [H A multi-domain (MD) cy- 
cloidal structure is described by averaging the calculated 
intensities for the CW and CCW spin structures. These 
results are compared to the measured polarized beam in- 
tensities in Fig. [3] for different domain population (SD 
and MD), field conditions (HF and VF), and spin- flipper 
configurations. Measured intensities in the absence of E 
agree with the calculated MD values, indicating equal 



CW and CCW volume fractions. The polarized beam 
data are in excellent agreement with the intensities cal- 
culated on the basis of the previously published structure 
0. The analysis accounts for the weak peaks in Fig. [2] 
and allows us to conclude that E applied along ±b in- 
duces a CCW(CW) cycloid, as defined in Fig. [l](b). 

We now examine the evolution of ferroelectric and cy- 
cloidal (i.e., multiferroic) domains with different electric 
field protocols, and show that these multiferroic domains 
in NVO have a defined magnetic and electric state that is 
directly coupled. Fig. [4] shows polarized diffraction and 
pyroelectric current data sensitive to the volume fraction 
of magnetic and ferroelectric domains, respectively, after 
E-cooling the sample. The magnetic diffraction data are 
presented in terms of the CW-CCW cycloid asymmetry 
{A): 

(/±) + (Jt) {a^^)-{a^cw)' 

Here I^{I±) is the magnetic Bragg peak intensity un- 
der the on/off (off /on) spin- flipper configurations and 
^cw^^ccw) calculated cross-section for CW 

(CCW) cycloids in the off/on configuration. This ex- 
pression correctly accounts for the contribution of each 
cycloid to /zp and I± . The data show that by cooling in a 
sufficiently high E oriented along the — b direction, 98 % 
of the long-range ordered part of the sample contains CW 
cycloids (the direction of E was confirmed in a separate 
experiment). The maximum spontaneous electric polar- 
ization attained was ~ 47 /iC/m^. This is 46 % lower 
than previously observed for NVO [5j]. The discrepancy 
may reflect uncertainty in determining the surface area 
for the small single crystalline sample or sample variabil- 
ity. 

Fig. [5] shows the E-dependence of the magnetic and 
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FIG. 5: (Color online) Hysteresis loop for magnetic asym- 
metry and electric polarization versus E. The sample was 
zero-field cooled to 5.0(2) K. Hysteresis data were obtained 
while varying E as indicated by the arrows. Typical waiting 
time before neutron diffraction measurements was 10 s and 
average ramp rate was 2 kV/(m-s). Ferroelectric hysteresis 
data were obtained with the sample in an evacuated can (to 
avoid electric breakdown in helium exchange gas) by integrat- 
ing the voltage-step induced current while varying E in steps 
of 1 kV/m at an average rate of 0.11 kV/(m-s). 



ferroelectric volume fractions after zero-field cooling the 
sample into the LTI phase. The product of coercive 
field and saturation polarization for NVO is 1.05 x 10~^ 
J-cm~^. This value is similar to multiferroic MnW04 
[13], but more than 4 orders of magnitude smaller than 
for the bulk ferroelectric Pb(Zn-Ti)03 and the thin- film 



ferroelectric BiFeOa [IJ, This reflects the differ- 



ent origin and much smaller magnitude of ferroelectric 
lattice distortion in inversion symmetry-breaking frus- 
trated magnets. Excellent correspondence between the 
E-dependent electric and magnetic domain fractions un- 
der field-cooled and zero-field cooled conditions leads us 
to conclude that these domains coincide in NVO. Fur- 
thermore, these results suggest that the magnetic and 
electric domain walls are coupled in this system. This 
is different from multiferroic hexagonal YMnOs, where 
magnetic and ferroelectric domain walls can move inde- 
pendently [IS]. NVO's multiferroic domains should have 
interesting domain wall characteristics that could be im- 
aged with magneto-optical techniques [17]. 

While the LTI cycloids could result from compet- 
ing Heisenberg interactions along the a axis, definite 
spin handedness throughout the sample as demonstrated 
here can only originate from an interaction like the 
Dzyaloshinskii-Moriya (DM) of the form D • (Sn x Sn+i) 
with D • a 7^ 0, which has opposite expectation values 
for the CW and CCW cycloids. The importance of DM 
interactions was previously pointed out for BiFeOs 0- 



Such a term is forbidden by symmetry in the paramag- 
netic phase of NVO [18|] , but apparently must be present 
in the Hamiltonian for the LTI phase. We infer that 
magneto-elastic distortions in this phase generate DM 
interactions that select a definite cycloidal handedness. 

In summary, we have demonstrated electric control of 
cycloidal handedness in multiferroic NVO through po- 
larized magnetic neutron diffraction. The previously de- 
termined spin structure is consistent with the polarized 
diffraction data and our quantitative analysis allows us to 
conclude that a single domain cycloidal structure can be 
generated through E-cooling to the level of 98 %. We 
observed comparable ferroelectric polarization and cy- 
cloidal magnetic hysteresis, indicating that electric and 
magnetic domains are coupled, as predicted by the trilin- 
ear theory. The results suggest that domains with elec- 
tric polarization along ±b feature DM interactions that 
select cycloids with CCW(CW) handedness. 
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